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When solid surfaces are immersed in aqueous solutions, some of their charges can dissociate and 
leave behind charge patches on the surface. Although the charges are distributed heterogeneously 
on the surface, most of the theoretical models treat them as homogeneous. For overall non-neutral 
surfaces, the assumption of surface charge homogeneity is rather reasonable, since the leading terms 
of two such interacting surfaces depend on the non-zero average charge. However, for overall neutral 
surfaces, the nature of the surface charge distribution is crucial in determining the inter-surface 
interaction. In the present work we study the interaction between two charged surfaces across an 
aqueous solution for several charge distributions. The analysis is preformed within the framework 
of the linearized Poisson-Boltzmann theory. For periodic charge distributions the interaction is 
found to be repulsive at small separations, unless the two surface distributions are completely out- 
of-phase with respect to each other. For quenched random charge distributions we find that due to 
the presence of the ionic solution in between the surfaces, the inter-surface repulsion dominates over 
the attraction in the linear regime of the Poisson-Boltzmann theory. The effect of quenched charge 
heterogeneity is found to be particularly substantial in the case of large charge domains. 



I. INTRODUCTION 

Long-range interactions between charged surfaces sub- 
stantially influence the structural properties of soft ma- 
terials, such as lipid membranes and colloidal suspen- 
sions [l|-|3( • The nature of the long-range interaction be- 
tween two charged surfaces immersed in an ionic solution 
is mainly determined by electrostatic interactions, medi- 
ated by the ionic solutes and polar solvent. When mod- 
eling the electrostatic interaction, surfaces are usually 
assumed to be homogeneously charged [l], How- 
ever, most charged surfaces in soft matter are heteroge- 
neous over a certain length scale. A schematic drawing 
of two such heterogeneously charged surfaces is presented 
in Fig. [TJ While the interaction between two homoge- 
neous surfaces depends on the average surface charge and 
the inter-surface separation, the interaction between two 
heterogeneous surfaces depends on the intra- and inter- 
surface charge correlations as well. 



Several experimental studies |7Hl7| measured inter- 
surface forces for different configurations of surface 
charge heterogeneities. For example, it was demon- 
strated 0-[IH that coating a negatively charged mica 
surface with a cationic surfactant monolayer may lead to 
the formation of positively charged bilayer patches, which 
neutralize the negative patches of the bare mica. A long- 
range attractive force was measured @-[H] between two 
such coated mica surfaces. In another setup, it was found 
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that adsorption of positive polyclcctrolytcs on negatively- 
charged spherical beads can lead to charge inversion of 
the beads 0, [I?} ■ The adsorbed polyelectrolyte forms 
positively charge domains, while patches of bare regions 
on the beads remain negatively charged. It has been 
shown [ijj that by controlling the amount of adsorbed 
polyelectrolyte, overall neutral particles with a patchy 
heterogeneous charge distribution can be prepared. In 
this special case of neutral beads, the inter-particle inter- 
action is attractive, while for smaller or larger amount of 
adsorbed polyelectrolyte (where the average bead charge 
is, respectively, negative or positive) the interaction is 
repulsive. 

The interaction between two periodically ordered 
charged surfaces has been studied theoretically in sev- 
eral works [l8l - |26| within the framework of the linearized 
Poisson-Boltzmann (PB) equation. For non-neutral sur- 
faces, the leading interaction term depends on the av- 
erage charge. For neutral surfaces, where the average 
charge is zero, the interaction depends strongly on the 
relative phase between the two charge distributions, and 
can vary from being repulsive to attractive. For two sur- 
faces with identical average charge and small amplitude 
charge modulation it was shown, in the non-linear regime, 
that the repulsive interaction is weaker than the repul- 
sion between two uniformly charged surfaces with the 
same average charge [271 32] . However, for non-zero av- 
erage charge, the charge modulation leads only to a small 
correction, as compared to the leading interaction term 
determined by the average charge. 

In a more general case of experimental relevance the 
surface charges are distributed randomly (and not pe- 
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FIG. 1: Schematic drawing of two planar heterogeneously 
charged surfaces immersed in ionic solution. The gray re- 
gions are positively charged and the white ones are negatively 
charged. The charge distributions of the bottom and top sur- 
faces are given by Eb (x, y) and Et {x,y), respectively. The 
inter-surface separation is d. 



riodically). We further make the important distinction 
between quenched and annealed cases of surface charge 
disorder. In the annealed case, the surface charges are in 
thermodynamic equilibrium with other system variables, 
such as inter-surface separation, temperature and con- 
centration of other species. The attraction can then be 
caused by self-adjusting of surface charge domains, where 
positively charge patches on one surface position them- 
selves against negatively charge patches on the second 
surface and vice versa [3314361 ] . For this case to hold, the 
typical time-scale of patch rearrangement must be sub- 
stantially shorter than the measurement time-scale, such 
that the self-adjustment (annealing) of charged domains 
on the surfaces has enough time to be completed. 

In the quenched case, the surface charges are frozen 
and independent of other system variables. For a 
quenched distribution of surface charges with no inter- 
and intra-surface correlations, it was reported that such 
heterogeneities have no effect on the two-surface interac- 
tion at the mean- field level [37], [HI , and the interaction 
depends solely on the average surface charge. 

Another interesting case of quenched disorder, to be 
considered in detail in the present work, is the case where 
finite-size charge domains (patches) are randomly dis- 
tributed on each surface. In some of the experiments 
such random charge domains stem from the specific sur- 
face preparation [7H15|. 

The outline of our paper is as follows. The model 
is formulated in Sec. [TTJ In Sec. IIIII the interaction be- 
tween surfaces with periodic charge modulations is re- 
visited and a crossover from attraction to repulsion is 
discussed. Quenched random charge distributions are 
treated in Sec. IIVI for several cases within the linear 
regime of the PB equation. In particular, the limit of 
small domains having a size in the molecular range is 
compared with the limit of large domains. In Sec. IVl we 
discuss the implications of our results, and concluding 
remarks arc presented in Sec. VI. Finally, in the limit of 
infinitely large charge domains, we compare in the Ap- 



pendix the predictions of the linear reg ime to those of the 
non-linear one considered in Ref. |15j . where an attrac- 
tive inter-surface interaction was predicted in the limit of 
infinitely-large charge domains. An interesting crossover 
is found from inter-surface attraction of strongly charged 
surfaces to repulsion for weakly charged surfaces, empha- 
sizing even further the role of the domain lateral size. 



II. THE LINEAR POISSON-BOLTZMANN 
MODEL 

We consider an aqueous solution as depicted in Fig. [TJ 
which is confined between two semi-infinite and planar 
charged surfaces, located at z = id/2, where the z-axis is 
perpendicular to the surfaces. The surface-charge density 
distribution (charge per unit area) is given by Eb (x, y) 
for the bottom surface at z = —d/2, and Et [x, y) for 
the top one at z = d/2. As we are interested in the 
effect of quenched charge heterogeneity, we model the 
lateral surface charge distributions Eb and Et as fixed 
charge boundary conditions. Namely, they are frozen and 
independent of other system variables. 

The aqueous solution contains a 1:1 monovalent salt 
ions. The solvent (water) is modeled as a homogeneous 
dielectric background with dielectric constant e = 80. 
The solution is coupled to a reservoir of ionic density 
nf,. The electrostatic potential, ip, and the ionic densi- 



ties, n± 



Teif'/k B T are ca i cu lated via the Poisson- 



Boltzmann (PB) equation: 



v 2 v> 



&irenb 



sinh (eip/ksT) , 



(1) 



where e is the electron charge, k^,T is the thermal energy, 
and /cb is the Boltzmann constant. By rescaling the po- 
tential ip — > 4> = eip/k^T the PB equation is rewritten 
for cj), the dimcnsionless electrostatic potential as: 



sinh « 



(2) 



where k d = A D = y/8irlBnb is the inverse Debye screen- 
ing length and Ib = e 2 / (efcsT) is the Bjerrum length. 
The boundary conditions at z = id/2 are given by: 



dz 



dz 



-d/2 



d/2 



(3) 



where a = AttIb (Ee/e) and rj = AttIb (St/e) are the two 
rcscalcd surface charge densities, having dimensions of 
inverse length. For homogeneous surface charge densi- 
ties the rescaled densities a and r\ are the inverse of the 
Gouy-Chapman length Q (up to a factor of two) that 
characterizes the thickness of the condensed ionic layer 
near a homogeneously charged surface. 
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The interaction energy between the two surfaces as a 
function of the separation d is given by the free energy, 
F(d), subject to the boundary conditions, Eq. (j3]). We 
obtain F (d) by using the charging method |J| , and cal- 
culate the work needed to increase the surface charge 
incrementally, at each point on the surface, from zero to 
the desired final value. The charging free-energy of two 
surfaces coupled with an ionic reservoir of density rib is 
given by: 



k B T 



F(d) 



r r f(<rfi) 

Kd // dxdy / da' ^[a'rf = G] Z= _ 
J J J (0,0) 



d/2 



+ n °JJ dxdy J o) dv ' ^^ ,T1 ^ z = d / 2 ' ^ 

where throughout the paper F (d) is made dimensionlcss 
by rescaling it, F — > (A-KK u l-Q/k-QT)F . The electrostatic 
potential on the surfaces, <f>, is written in the above equa- 
tion as a function of the distributions a and rj to stress 
that the integration is performed with respect to these 
variables. 

Equations (l)-(4) are valid for the general PB theory, 
while hereafter we concentrate on the linear PB equation, 
being valid as long as = eip/k-oT <C 1 (ip -C 25mV): 



V 2 = ^0. 



(5) 



The in-plane Fourier transform of <f>(p,z) in the p = 
(x, y) plane is denoted by </>k (z) and is given by: 



0k(z)= / d>e ik ^(p,z) 



(6) 



where k is the inplane wavevector. Using k we get from 
Eq. ©: 



dz 2 



(7) 



The boundary conditions for k are similarly obtained 
by Fourier transforming Eq. 



0z 



-d/2 



0z 



?7k 



(8) 



d/2 



where a^, r/k are, respectively, the Fourier transforms of 
a{p) and rj(p): 



a k = J d 2 p^»a{p) , 



Vk = 



d 2 pe iU - p T 1 (p) . 



(9) 



Solving for the z-dcpendcncc of 0k, Eqs. (O-© yields: 
(>k + ?7k) cosh (qz) (r/ k - cr k ) sinh (qz) 



4>k (z) 



2q smb. (qd/2) 



2q cosh (qd/2) 



(10) 

The free energy can be re-expressed in terms of the 
Fourier transforms <7k, r) k and 0k- By substituting 
Eq. (fTU)) into the free energy, Eq. Q, we obtain: 



F(d) = 



d 2 k <7 k ?7-k + cr_k?/k + (<7k0-_ k + VkV-k) c 



— qd 



2q sinh (qd) 



(11) 



where the reference contribution at d — s- 00 was sub- 
tracted in order to obtain zero interaction between sur- 
faces at infinite separation, F(d — > 00) = 0. Because of 
the linearity, the PB free energy can also be written as a 
sum of two decoupled variables = (<7k ± 77k) 1 2: 



-qd/2 



q sinh (qd/2) 

Q -qd/2 

gcosh (qd/2) 



'k"-k 



(12) 



The first term is positive definite and decreases monoton- 
ically as a function of the separation d, contributing to 
the repulsive part of the inter-surface interaction, while 
the second term, being negative definite and a monotonic 
increasing function of d, gives the attractive part of the 
interaction. 

Quite generally, it can be stated that the magnitude 
of the repulsive term is larger than the attractive one 
for all values of q and d, because of the inequality 
sinh (qd/2) < cosh (qd/2). Furthermore, for d — !> the 
magnitude of the repulsive term diverges, while the at- 
tractive term approaches a finite limiting value. Conse- 
quently, only when cr^ k = for any wavenumber k, 
does the interaction become attractive for any d. This 
occurs only when the distributions of the two surfaces 
are completely out-of-phasc (anti-symmetric) with each 
other, a (p) = —r](p). Furthermore, by properly choos- 
ing cr^, F(d) can be made attractive for large enough 
d. 



III. PERIODIC SURFACE CHARGE 
MODULATION 

In the previous section the distributions a (p) and rj (p) 
have been taken to be arbitrary. Hence, the inter-surface 
interaction of Eqs. ([TTj) - (|12p applies to any form of sur- 
face charge distributions. We treat next two interacting 
charged surfaces each characterized by a periodic charge 
modulation. As periodic charge modulations have been 
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investigated previously [18H26I ]. we revisit the periodic 
case in order to discuss the interesting crossover from at- 
traction to repulsion and the limiting behavior for small 
and large periodic domains. It is also instructive to com- 
pare the periodic case with the quenched disorder that 
will be derived later in Sec. IV. 



A. Single Mode 

We first consider two surfaces characterized by the 
same singlc-kg mode modulation 



cr (p) = C a cos (k • p) , 



77 (p) = C v cos (k • p + ip) 



(13) 



where C a and C v (having units of inverse length) are 
the modulation amplitudes of the two surfaces, taken 
to be positive without loss of generality. The angle 
— 7r < if < 7r is the relative phase between the two charge 
modulations. It reflects an in-phase arrangement of like 
charges at ip = 0, and an anti-phase one at ip = ±tt. 

The inter-surface interaction per unit area, F (d) /A, 
is given by: 



F (d) _ 2k u C C^ cos ip + K D c-«> d {Cl + C*) 



A 



4q sinh(god) 



(14) 



k Q . 



where A is the lateral surface area and qo = 
For small separations, qod <C 1, the energy scales as 
F I A ~ « D (<?o^) 1 while for large separations, qod 3> 1, 
the energy decays exponentially, F/A ~ n D qQ 1 e~ g ° d . 

Rewriting F/A in terms of C± = (C a ±C n ) /2 as in 
Eq. dHJ leads to: 



F(d) 
A 



-q d 



2q sinh (q d) 

K 



Ci 



e -q d _ cog ^ 



2q sinh (q d) 



Ci 



(15) 



Here both the C\ and Ci terms can be either repulsive 
or attractive, depending on the values of d and p. 

In Fig. [2] the inter-surface interaction per unit area 
F (d) I A is presented for charge modulations with two 
equal amplitudes, C G = C n = C+ . The second term of 
Eq. (jT5"j) then vanishes as C_ = 0, and the sign of the 
interaction energy is determined by the relative phase ip 
and d. The repulsion is maximal with respect to tp for 
the in-phase state p = 0, while for the anti-phase state, 
ip = ±7r, the attraction is maximal: 



F(d) 
A 



2qo sinh(go^/2) "+- 



2q a cosh(q d/2) ° + 



for ip = 



for p = ±7r 



(16) 



Furthermore, for modulated surfaces with equal ampli- 
tudes, C a = C^, the interaction in Eq. (|15|) is purely 
repulsive in the range \tp\ < 7r/2, where the cosine is 
positive. When cos ip becomes negative, in the range 
7r/2 < \ip\ < 7r, the inter-surface interaction F (d) has 
a minimum at d c = q^ 1 coshT 1 (1/ \cosp\), leading to 
a crossover from attraction for d > d c to repulsion for 
d < d c . 

The interaction between two non-neutral and uni- 
formly charged surfaces is obtained by setting ko = p = 
in Eqs. (fT2| - ([T3|) . recovering the well-known result of 
Parsegian and Gingell [39l |: 



< D d/2 



F(d) 
A sinh(K D d/2) 



-K d/2 
e q2 

cosh (n D d/2) 



(17) 



written here in a decoupled way, separating the pure re- 
pulsive from the pure attractive contributions in terms of 
the amplitudes C± — (C a ± C v ) /2, where C a and C v are 
uniform surface charge densities. Unlike neutral modu- 
lated surfaces, for two uniformly charged surfaces with 
the same sign (C' a ■ C v > 0), the interaction is purely 
repulsive since the inequality C + > C_ is always valid. 

For uniform charged surfaces, only the ions in solution 
screen the interaction, whereas for modulated charged 
surfaces there is an additional screening effect, yielding 



an effective screening length, q 1 = 1/y k'^ + k^. There- 
fore, the lateral charge modulations result in a faster de- 
cay of the interaction as compared to the interaction be- 
tween two uniformly charged surfaces, due to the addi- 
tional mechanism of screening. 

Note that when two surfaces have each a different 
single- mode modulation, k a 7^ k^, the integral over the 
coupling terms in Eq. (fTTj) vanishes, leading to a purely 
repulsive interaction: 



-q<jd 



F(d) _ * D 
A Aq a smh.(q a d) 



C 2 



jr, sinh(q v d) v ' 



(18) 



are 



where q a = ^/k^ + fc 2 ., q v = J + k%, and k a , k v 

the wavenumbers of the bottom and top surfaces, respec- 
tively. 

The inter-surface interaction as expressed in Eq. (fT5)) 
can be interpreted as a superposition of two decoupled 
systems each composed of one uniform charged surface 
and a second (virtual) surface that is neutral. Namely, 
the first term corresponds to a uniformly charged sur- 
face located at z = —d/2 and a neutral one located at 
z = +d/2, while the second term corresponds to a neutral 
surface located at z — ~d/2 and a uniformly charged sur- 
face located at z = +d/2. Although the virtual surface is 
neutral and docs not interact electrostatically, it leads to 
a repulsive interaction because the counter-ions are con- 
fined between the two surfaces and their translational 
entropy is reduced. The magnitude of the confined ionic 
volume depends linearly on d, and decreases at smaller 
separations. 
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FIG. 2: The inter-surface interaction (per unit area) F/A for a single go-mode, rescaled by the prefactor of Eq. (|15[) , F — 
(qo/K D C+) F, and plotted for the equal amplitude case: C a = C v = C+ (C_ = 0). (a) The inter-surface interaction as a 
function of the inter-surface separation d in units of q^ 1 . The values of the relative angle <p are 0, 0.85-7T, iv for the solid, dashed 
and dotted lines, respectively, (b) The inter-surface interaction as a function of the relative phase (p. The values of qod are 0.5, 
1 and 2 for the solid, dashed and dotted lines, respectively. 



The above result suggests that a slight difference be- 
tween the two wavelengths (k a ^ k n ) is sufficient to elim- 
inate any inter-surface attraction between two periodi- 
cally ordered surfaces, implying that the attraction is a 
delicate effect that might be difficult to obtain in experi- 
mental conditions corresponding to the linear PB regime. 



B. Multi-mode Modulation 

The analysis of the previous section can be generalized 
to multi-mode distributions of the charges on the two 
surfaces. The charge distributions are written as: 



6a (p) 
5r] (p) - 



■ 00 

Vo, 



(19) 



where <ro = ( cr ) an d ?7o = (v) are t ne average surface 
charge densities, and 6a (p) and 6r] (p) are the charge 
modulations around these averages. The inter-surface 
interaction of Eq. (jTTJ) can be separated into two terms 



F (d) =F (d) + F 1 (d) , 



(20) 



where the first term depends only on ao and 770, and 
coincides with Eq. (fT7)) with C± = (ao ±r]o)/2: 



F (d) = 2a r, + (<xg + ?fa 2 ) C -^ d 
A 2sinh(K D d) 



(21) 



It vanishes when both surfaces are neutral, a = t}q = 
0. The second term is the contribution due to charge 



heterogeneity, given by Eq. (fTTj) without the mode k = 0: 



Fi (d) = K D 



d 2 k 



Sa k 6r]_ k + 6a-u5r)k 
k#o (2tt) 2 L 2q sink (qd) 

(da k 5a- k + 6ii k 5i]_ k ) c" 
2q sinh (qd) 



(22) 



where 5a^ and Sr] k are the Fourier transforms of da (p) 
and Sr] (p) , respectively. By definition, from Eq. (fTT)]) , 
(6a) = (6rf) = 0, or equivalently in Fourier space <5<7k=o = 
< ^ 7 ?k=o = 0- Note that Sa^ and 6r] k have units of length 
in Fourier space (unlike the units of C a and C v ). Hence, 
F\ is dimcnsionless. 

An important length scale for a general periodic 
(multi-mode) distribution is periodic charge domain size, 
L, which is related to the smallest wavenumbcr, L = 
27r/fc m i n . The wavevector k m ; n also acts as a lower cutoff 
in the k-spacc integration, Eq. (|2"2"1) . 

Two limiting cases of small and large fc m i n can be con- 
sidered separately. In the limit of charges with molecular- 
size heterogeneities, the domain size is molecular L <C 

k~\ and <j min = ^k^+ fc^ in ~ 2ir/L > k d depends 

inversely on L. Therefore, the leading term in the inte- 



grand varies a,s n D q c 



k Lc 27Td / L (as long 



as L <§; d), and leads to a negligible contribution due to 
charge heterogeneity, F\ (d) — > 0. 

On the other hand, in the limit of large domains, L 3> 
K" 1 , g m i n ~ n D + O (L~ 2 ) has no dependence on L. The 
decay length of the interaction is the Debye length, k" 1 , 
similarly to the average charge contribution in Eq. (|21[) . 
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The leading term of Fi(d) is then given by: 



Fi(d) 



Si + S 2 e" 



sinh (n D d) 



(23) 



where the prefactors 

S i = 5 / (f^j2 (^k^-k + dv-kSrik) , and S 2 = 

1 1 (trr) 2 (l^ g k| 2 + |<S?7k| 2 ) > do not depend on d, and 
can be interpreted as contributions from effective uniform 
surface charges. 

The similar dependence on d of Eqs. ([21]) and ([23| im- 
plies that the effect of large domains can be understood 
as a modification of the prefactors in Eq. (|2"Tj) . When Si 
is negative, the contribution F\ (d) ~ 2Sie~ K o d < 0, is 
attractive for large separations K^(i> 1, and is followed 
by a crossover at d c = k" 1 cosh" (— S2/S1) from attrac- 
tion for d > d c to repulsion for d < d c , resembling the 
crossover found in Eq. (| 15[) . 

Examining the two limits of small and large domains 
leads to the observation that periodic charge modula- 
tion has a significant effect only when the domain size 
is sufficiently large. Thus, for charge modulations char- 
acterized by molecular heterogeneities, one can treat the 
interaction effectively as if the charges on the two sur- 
faces are distributed homogeneously. On the other hand, 
for surfaces where the charge domains are larger than 
the Debye screening length k^ 1 , the contribution due to 
charge modulation cannot be neglected, and may lead to 
a substantial change in the strength of the inter-surface 
interaction. Furthermore, a crossover from repulsion at 
d < d c to attraction at larger separations, d > d c , can be 
induced. 



Time-dependent Lateral Displacement and 
Inter-surface Correlations 



The case of two identical single-mode distributions 
with a relative phase (p, Eqs. (jT3 |) - fp~5]) . can be gener- 
alized for two identical multi-mode distributions with a 
relative lateral displacement £: 



6r) (p + £) = da (p) 



(24) 



The above charge distribution, Eq. (|24|) . is motivated by 
a specific inter-surface force experiment setup [l5[ that 
will be discussed in Sec. V, while here we explore its 
theoretical consequences. The relative displacement £ = 
(&x,£y) is related to a fc-dependent phase ipt- = k • £ in 
Fourier space with <5r?k 
reads: 



5<r k c lipk , and F 1 from Eq. <J22 



Fi (d) = « D 



d 2 k cos (k ■ £) + e~i d 
(2?r) 2 gsinh (qd) 



(25) 



For £ = the two distributions are in-phase (symmetric) , 
leading to repulsive contributions from all modes for any 
d. In the more general case of a relative displacement 
£ 7^ 0, each k-mode can have either an attractive or 



a repulsive contribution, and Eq. (|25|) gives the overall 
positive/negative sign of Fi, due to strong correlations 
between the two surfaces. 

We consider two surfaces with identical charge dis- 
tributions that are relatively displaced by applying on 
one of them a time-dependent lateral force. This force 
leads to a time-dependent lateral vector £(t) such that 
5r](p + £(t)) = Sa(p). The time dependence £(t) is re- 
lated to the temporal dependence of the applied lateral 
(shear) force. The validity of this model holds as long 
as the lateral motion is a quasi-static process, for which 
the ions in the solution equilibrate faster than the typical 
time-scale of the surface lateral motion. Then, the inter- 
surface interaction is given by averaging the cos [k ■ £(t)] 
term in Eq. (|25[) over the displacement period T : 



(Fi (d)) 5 



d 2 k (cos [k -£(t)] ) T + e 



-qd 



(27T) 



q sinh (qd) 



(26) 

where F\ is the contribution to the inter-surface in- 
teraction due to charge modulation, and (0) T = 

T _1 J Q T dtO(t) is the time average. Since the average 
of the cosine depends on the lateral periodic motion, the 
interaction energy would vary when changing the oscilla- 
tory mode of the force. 

A simple example is a square-wave motion in the x 
direction, given by: 



£(0 < t < T) 



(£ - A£) x 
(£ + M) x 



< t < T/2, 
T/2<t< T, 



(27) 



where £q is the mean lateral displacement between the 
two surfaces, and A£ is the oscillation amplitude. The 
inter-surface interaction is given by 



(F 1 (d)) T = 



d 2 k cos (k x £ ) cos (k x A£) + e~ qd 



q sinh (qd) 



l&Tkl 



(28) 



and depends both on the mean displacement £ and the 
oscillation amplitude A£. Here k x is the wavevector com- 
ponent along the direction of the shear motion. 

For a general periodic dependence £(t), if the ampli- 
tude of the lateral oscillations A£(t) = £(t) — £q is small, 
k x A£ <C 1, and to leading orders the interaction depends 
on the average displacement, £q, and the mean square 
amplitude of the oscillation, (A£ 2 )t = A£ 2 : 



(Fi (d)) q 



d 2 k 
(2^) 2 



(l-ifc 2 (A£) 2 )cos(Mo) 



q s'mh(qd) 



q sinh (qd) 



|^k| 



(29) 



Note that the above equation is exact for the square- wave 
distribution of Eq. (J27J. 
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FIG. 3: The averaged interaction energy, F/A as a function of the inter-surface separation d for randomly charged surfaces. 
In (a) the free energy from Eq. ((33}, for an uncorrelated Gaussian distribution Eq. (|32[) . is plotted on a semi-log plot as a 
function of n D d. In (b) the free energy is plotted on a semi-log plot as a function of «; D d, where the charge distribution is 
a two-dimensional Lorentzian, Eq. (|36[). Three values of k d £ are shown: k d £ = 100 (solid line), k d £ = 1 (dashed line), and 
k d £ = 0.1 (dotted line). In the inset the dependence of the free energy on k d £ is plotted on a log-log scale for two values of the 
separation, n D d — 0.1 (solid line), and K D d = 1 (dashed line). The plotted free energy is rescaled by the prefactor of Eq. (|33p 
F = 2tt (F) I (n%'fa 3 ) for (a), and F = (F) / 7 2 from Eq. (36]) for (b). 



IV. QUENCHED SURFACE DISORDER: 
PATCHY SURFACES 



We generalize now the periodic distribution of Sec. Ill 
to patchy surfaces with random charge domains, but 
which are overall neutral. Consider two surfaces with 
charge distributions a (p) and rj(p), respectively. The 
surface charges arc randomly distributed with a joint 
probability distribution P [a (p) , r\ (p)]. The inter- 
surface interaction, Eq. is obtained by averaging 
over the bilinear terms of the Fourier transform: 



d 2 k 



(7,7} 



<T,T) 



<T,T) 



2q sinh (qd) 



-qd 



2q sinh [qd) 



(30) 



where (0) ar] = J Va J VrjP [a,r]]0 is the average over 
the joint a and i] distribution. In experiments the prob- 
ability distribution is usually determined by the prepa- 
ration procedure of the surfaces. 

When the surfaces are prepared separately, there are 
no inter-surface correlations, (ck^-k) = 0, leading to a 
purely repulsive interaction: 



d 2 k (o- k o-_k) 



2q sinh (qd) 



'_v e - qd 



>0, (31) 



and the strength of the repulsive interaction depends on 
the probability distribution of each surface. For an uncor- 
related Gaussian distribution, the two-point correlation 



function is given by: 

(a(p)a(p')) a = { V (p) V (p% = 7 2 5(^^\ , (32) 

where S(p) is the two-dimensional Dirac (5-function, 7 
is the root mean square charge density (taken to be the 
same on the two surfaces), and a is a conveniently defined 
molecular length. This leads to (crkcr^k)^ = (VkV-k) v = 
A^a 2 , and to an inter-surface interaction per unit area: 



(F(d)) c 



K ^ 2 a 2 



where 



I( Ku d) 



dk 



2- 



ke~i d 



I{n D d) > , 



ln(l 



(33) 



n D q sinh (qd) 



K„d 



(34) 

The repulsive free-energy for the uncorrelated Gaussian 
case is presented in Fig. [3fa) . At separations larger than 
the Debye screening length, k d ((> 1, the leading term 
decays exponentially as (F) /A ~ c~ 2fi o d / (n D d), while 
at small separations, K v d <C 1, the energy scales as 
(F) /A~In[l/(2« D d)]/(« D d). 

When calculating these scaling relations we ignored the 
molecular length scale a, since I(K D d) has no dependence 
on a. However, the prefactor in Eq. (f3"3"|) scales as the 
square of the molecular length scale (F)/A ~ (n D a) 2 . 
Hence, the contribution of uncorrelated Gaussian dis- 
order, Eq. (|32p . vanishes in the limit where a is much 
smaller than the Debye length, K n a <C 1, in agreement 
with the results reported in Ref. [371 ]. 

In order to model the possibility of finite (macroscopic) 
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charge domains we replace Eq. (|32j) by a two-point cor- 
relation function having a Lorentzian distribution: 



i 2 e 



(p-p'Y+e 



(35) 



where £ is the charge correlation length on the sur- 
face, which can be associated with the domain size. In 
Fourier space this leads to (crkf-k)^ = (??k?7-k) n = 
2nA'y 2 ^ 2 K (fc£), where K (x) is the zeroth-order modi- 
fied Bessel function of the second kind. 

The inter-surface interaction is then given by: 



( F W)*,r, 2,2 



.4 



dfc^MMl, (36 ) 
q sinh (qd) 



and the dependence of (F(d))„- tri on n D d for three val- 
ues of the correlation length £ is shown in Fig. [3jb). In 
the limit of short-range correlations [dotted line in[3]Jb)], 
k d £ <C 1, the integrand in Eq. (|3T)|) depends weakly 
on £, and the interaction energy scales as (F) j A ~ 
— ln(K D £)£ 2 , resembling the scaling of the Gaussian dis- 
tribution case, Eq. (|3"3")h where (F) ~ a 2 . In the other 
limit of large domains [solid line in|3{b)], where k d £ 1, 
the main contribution to the integral comes from the 
small values of fc, and the interaction is independent of £, 
(F) j A ~ or K v d j sinh (n D d). This reproduces repulsion 
with a decay length that does not depend on the lateral 
length scale, £. It is effectively analogous to repulsion 
between uniform surfaces as in Eq. (|17p with an effective 
separation that is twice as large as that of the homoge- 
neous case. These two limiting behaviors of (F) /A as a 
function of k d £ are shown in the inset of Fig. ^b) for 
two values of d. 

By examining the short- and long-range charge distri- 
butions we find that short-range charge disorder has a 
negligible effect on the inter-surface energy, while long- 
range disorder leads to a substantial effect, which cannot 
be omitted even in the case of non-zero average charge. 
This observation is in accord with the results presented 
in the multi-mode periodic case, Sec. IIIIB1 where it was 
found that charge distributions with large charge do- 
mains changed the inter-surface interaction substantially, 
while distributions with small domains lead to a negligi- 
ble effect. 



V. DISCUSSION 

The model we present accounts for the inter-surface 
interaction of two heterogeneously charged surfaces. It is 
formulated on the mean-field level where the electrostatic 
potential is given by the Poisson-Boltzmann (PB) equa- 
tion. Our main aim is to study the effect of quenched 
charge disorder with finite domains on the inter-surface 
interaction between the two surfaces. 

We focus on the linear regime of the PB equation, 
where it is found that the inter-surface interaction of two 



arbitrary charge distributions depends on a bilinear cou- 
pling between the Fourier components of the two distri- 
butions. It is shown that the interaction can be expressed 
as a sum of an attractive term and a repulsive one. The 
repulsive term at small separations is usually larger than 
the attractive one, leading to a repulsion for d — > 0. In a 
special case where the two distributions are exactly out- 
of-phasc, a (p) = — ?/(p), such that negatively charge 
domains are located against positively charge domains in 
the other surface and vice versa, the interaction is attrac- 
tive for any value of d. 

Several types of surface charge distributions are stud- 
ied in the linear PB regime. First, we consider two over- 
all neutral surfaces with the same single-mode periodic 
charge density, but with a relative phase shift tp between 
the two. When the average surface charge is taken to 
be zero, the inter-surface interaction stems only from the 
lateral charge modulation. When the top surface mode 
is identical with the bottom one, the interaction can vary 
between pure attraction and pure repulsion as function 
of the relative phase ip between the two surfaces. For 
an out-of-phase distribution, ip = w, there is a relative 
displacement by half a wavelength between the bottom 
and top surfaces, and the interaction is purely attractive, 
while for the in-phase case, ip = 0, it is purely repulsive. 

These two limits are a consequence of strong inter- 
surface correlations. For intermediate values of the 
relative phase <p>, the correlation is weaker and a 
crossover from attraction (at large separations) to re- 
pulsion (at small separations) occurs at a separation 
d c ~ cosh^ 1 (l/|cosy|). Note that the crossover sepa- 
ration d c does not depend on the modulation amplitude. 
In the more general case where the bottom and top sur- 
faces have different (single) fc-mode, the interaction is 
found to be purely repulsive. The origin of the repulsion 
in this case is due to the confinement of the counter-ions 
in between the surfaces. 

The role of the domain size in a general periodic (multi- 
mode) distribution is also investigated. It is found that 
the contribution due to periodic charge modulation is 
important only in the limit of large domains k d L 3> 1. 
For sufficiently large domains the contribution due to 
charge heterogeneity is substantial, even in the case of 
non-neutral surfaces, where the main contribution stems 
from the average surface charge. 

The time-dependent relative displacement case dis- 
cussed in Sec. III.C is motivated by recent surface force 
experiment [l5[ , which examined whether the charge het- 
erogeneity on the two surfaces is annealed or randomly 
distributed (quenched case). In the experiments, the nor- 
mal forces between two mica surfaces partially coated 
with a cationic lipid bilayer were measured, while during 
the vertical approach the surfaces were also sheared later- 
ally in an oscillatory mode at a rate which is slower than 
the typical time required for the ions in solution to rear- 
range, but faster than the time scale of vertical approach 
or lipids rearrangement on the coated surfaces. Thus, it 
is assumed that the initial ordering of the surface charge 
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on each surface is preserved during the normal approach. 

If the hypothesis that the surfaces are periodically or- 
dered [34| is correct, then the measured inter-surface 
force should depend on the average and mean square am- 
plitude of the lateral motion as in Eqs. ([28|) - ([29| . How- 
ever, it was shown in Ref. |l5[ that the forces measured 
for different lateral motions remained unchanged, sug- 
gesting that the surface charge distribution is not peri- 
odic, but presumably is quenched and random. 

While in the annealed case the surface charges can re- 
arrange themselves, and the system approaches the min- 
imum energy state as was investigated in Refs. [33T - l36j ] . 
in the quenched case the accessible configurations of the 
surfaces are frozen. The quenched probability distribu- 
tion is usually determined by the experimental setup. 
For example, in surface force experiments the two sur- 
faces are prepared separately before the force measure- 
ment, and their corresponding charge distributions can 
be assumed to be independent with no inter-surface cor- 
relation. We consider this case and find that the in- 
teraction is purely repulsive in the linear PB model, 
due to the vanishing of inter-surface bilinear coupling 
terms. The inter-surface interaction can be thought of 
as an average over equal weights of repulsive and at- 
tractive contributions. However, since the repulsive term 
varies as exp(—qd/2)/ smh(qd/2) and the attractive one 
as exp(— qd/2)/ cosh(qd/2), the repulsive contribution is 
always larger than the attractive one, leading to an over- 
all repulsive interaction between uncorrelated randomly 
charged surfaces in the linear PB model. 

We discuss two specific types of random distributions. 
For Gaussian charge distributions with no inter- and 
intra-surface correlations, the interaction depends on a 
molecular length scale a (lower cutoff length). Disorder 
at the molecular level may emerge due to local processes 
that are not affected by neighboring surface charges. The 
contribution of the heterogeneity in this case is negligible, 
in agreement with Ref. [37] . 

In a more general case of Lorentzian distribution with 
correlation length £, we find two limiting regimes. For 
large values of £ we obtain a substantial effect due to 
charge disorder, while for small values of £ the limit of 
disorder at the molecular level is recovered and the con- 
tribution due to disorder is negligible. Note the similarity 
between these two limiting regimes and the case of peri- 
odic distributions mentioned earlier. For the latter and 
in the limit of large domains, nuL 3> 1, the effect is sig- 
nificant, while in the small domain limit, k^L <C 1, the 
domain-size effect is quite negligible. 

All the results obtained in the present work assumed 
a linear PB regime and demonstrate that quenched 
charge disorder leads (except in special set-ups) to a non- 
vanishing repulsive interaction between charged surfaces. 
It is reasonable to expect that similar considerations can 
be extended to the non-linear PB regime, where the elec- 
trostatic potential is too large for the linear approxima- 
tion to be valid. However, using a simplified model for 
the limit of infinitely-large charge domains, it was shown 



in Ref. [15( that the inter-surface interaction in the non- 
linear regime may become overall attractive. 

In the Appendix we reproduce this simplified treat- 
ment and compare its results, for a range of surface 
charge densities, with those of the linear model presented 
above. For weakly charged surfaces, when the validity of 
the linear PB approximation can be justified, the over- 
all repulsion predicted in the linear regime agrees with 
the non-linear one, while for strongly charged surfaces 
the linear PB equation fails to predict the attraction ob- 
tained in the non-linear regime. This suggests a possible 
crossover from repulsive to attractive inter-surface inter- 
action, as a function of charge density and patch size. 



VI. CONCLUDING REMARKS 

The results presented in our work demonstrate that 
charge heterogeneity may have important implications on 
the inter-surface interaction both for periodically modu- 
lated surfaces and randomly quenched ones. Several ex- 
tensions to our work can be considered. 

An effect not taken into account in the present work 
and left for future investigations is to look more carefully 
at the structure and thickness of the charged surface. 
This may be of importance for charged membranes com- 
posed of mixtures of charged and neutral lipids. As the 
lateral membrane charge density may be correlated with 
structural undulations and variation of membrane thick- 
ness, the inter-surface interaction between two mem- 
branes might be affected. 

A further remark should be made on the divergence of 
the free energy at small separations. This is related to 
our simplifying assumption on the boundary conditions, 
Eq. <j3j> , in which the electrostatic field is confined in be- 
tween the surfaces, and does not leak to the outer region. 
Namely, the electrostatic field in the outward direction of 
the surface is set to zero. This assumption is reasonable 
when the dielectric constant of the outer region is much 
lower than the one of the ionic solution. However, for suf- 
ficiently small inter-surface separations, one might need 
to consider the leakage of the field to the outer region 

Since the attraction in the non-linear regime was pre- 
dicted only in the limit of infinitely-large surface patches 
fl5j , it is worthwhile to investigate further the non- linear 
PB model. It will be of interest to study the origin of 
the attraction and relate it to the difference between the 
scaling of the counter-ion entropy for counter-ion release 
(attraction) versus counter- ion confinement (repulsion). 

Another possible extension would be to find a suitable 
approximation for the non-linear PB regime in three di- 
mensions with heterogeneous boundary conditions. Us- 
ing such an approximation may shed more light on the in- 
teraction between randomly charged surfaces within the 
PB theory, and related experiments. 

Acknowledgements. Wc thank Henri Orland and 
Rudi Podgornik for useful discussions and suggestions, 
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FIG. 4: The patch-patch inter-surface interaction energy per unit area, |-F|/A as a function of the separation d rescaled by 
k d . The plotted free energy F = F/k^ is rescaled by k^. (a) Attraction |-F a tt| (dashed line) and repulsion F Icp (solid line) 
are calculated by the non-linear PB equation for a — —rj = 10k d and a = r\ — 10k d , respectively. The inset shows the 
average of the repulsive and attractive corresponding free energies F av = (Frcp — |fatt|)/2, calculated by the linear and non- 
linear PB equation, (b) The solid and dashed lines are calculated with the non-linear PB equation for a = rj = 0.5k d , and 
a — —rj = 0.5k d , respectively. The inset shows the average F av . The results of the linear PB equation are omitted in (b) since 
their difference from the non-linear calculation is invisible. 
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Appendix: Comparison of the Linear and Non-linear 
PB Theory for Infinitely-large Charge Domains 

We repeat here the calculation presented in Ref. [l5[ for 
the non-linear PB theory and compare it with the linear 
theory treated in our paper. We consider two heteroge- 
neous surfaces where the typical size of a charge patch 
(domain) is much larger than any other length scale in 
the system. Each patch can be treated as an extended 
section of a uniform charged surface, and the leading in- 
teraction term would presumably be a superposition of 
all patch-patch interactions between the two surfaces. 

The two surfaces are overall neutral, and on each of 
them there are positive and negative charge patches, with 
uniform charge densities: a = ij = ±<7 , leading to four 
types of possible patch-patch interactions between the 
two surfaces, (±00 7 ^ za o)- The electrostatic potential ip is 
obtained by solving numerically the non-linear PB equa- 
tion for two uniformly charged surfaces. The patch-patch 
interaction energy is repulsive for the like-charge case, 
-Frep > 0, and attractive for the opposite-charge case, 
-Fatt < 0, and is obtained by integrating over the osmotic 
pressure H(d), 



/>OC 

F/A = 4:TTK D l B /k B T dzU 

Jd 



where the osmotic pressure is given by 
2 



8tt V dz 



-2k B Tn b 



cosh 



eif) 



1 



(A.2) 



(AT) 



and has the same dependence on ip for both attrac- 
tive and repulsive cases. The overall free energy F av 
is then given by averaging over the four possible patch- 
patch arrangements, yielding F av = (F lep + F att ) /2 = 
(F rep -|F att |)/2. 

In Fig. |4ta) a comparison between the free energy of 
the two cases is presented for strongly charged patches, 
a = ±J] = 10k d . At separations larger than K D d > 0.043, 
the attractive interaction (oppositely charged surfaces, 
(7 = — 77, dashed line) is stronger than the repulsive one 
(equally charged surfaces, a = r), solid line). This effect 
is related to different scaling of the ionic entropy as a 
function of d for the repulsive and attractive cases. At 
smaller separations, K v d < 0.043, the repulsive interac- 
tion becomes stronger than the attractive one, and even 
diverges. The inset shows the overall inter-surface inter- 
action as calculated in the simplified non-linear and linear 
PB models. It is clear that for large values of a and rj 
the linear approximation fails to predict the crossover be- 
tween attraction to repulsion of the non-linear PB equa- 
tion. 

The origin of the difference between the linear and non- 
liner PB models stems from the different dependence of 
the counter-ions entropy on the inter-surface separation d 
in the case of attraction and repulsion @, Ho| ■ While the 
entropy in the repulsive case is determined by counter-ion 
confinement 0] , the entropy in the attractive case is gov- 
erned by counter-ion release (40| . By directly comparing 
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the linear and non-linear PB solution for strongly and 
equally-charged surfaces, it is found that the counter-ion 
entropy is largely overestimated in the linear PB solu- 
tion [Fjj . This is due to a slower decay of the ionic profile 
close to charged surfaces in the linear case. 

For complementarity, in Fig.|4jb), we present the inter- 
surface interaction of weakly charged surfaces a = ±77 = 



0.5k d . As expected, the linear approximation agrees well 
with the simplified non- linear PB, and the overall interac- 
tion is repulsive. These rough estimates suggest a possi- 
ble crossover as a function of patch charge strength from 
repulsive to attractive inter-surface interaction, but re- 
quires further investigations of the non-linear PB equa- 
tions for surfaces with finite-size charge patches. 
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